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Abstract 

C++ package BtoVVana is aimed to extract physics parameters in decays 
of the scalar neutral B-meson to two vector particles in the intermediate state, 
and two leptons and two pseudoscalar mesons in the final state. Improved an- 
gular moments method is implemented in the package BtoVVana: the method 
uses the time dependent weighting functions and time-integrated obscrvables 
with variable upper time limit. These observables allow to keep full time in- 
formative contents of the decay. By using the package BtoVVana it was shown 
that the statistical errors of extraction of the observables strongly depends on 
the choice of the angular weighting functions. The best angular weighting func- 
tions, which are the linear combinations of the amplitude angular functions, 
are implemented in the package BtoVVana. The first version of the program 
contains algorithms to extract physics parameters in the case of untagged de- 
cays. The package BtoVVana and program codes for the channels J /'4'4' 
and J /iJjK* in the Monte Carlo generator SIMUB were written in parallel 
and mutually tested with high precision. 







Program summary 



Title of the program: BtoVVana 
Catalogue identifier: XXXX 
Program summary URL : XXXX 

Program obtainable from: http:/ /cmsdoc.cern.ch/~ shulga/BtoVVana/Bto VVana.html 
Computer: PC, two Intel 2.0 GHz processors, 512MB RAM 

Operating system: Linux Red Hat 6.1, 7.2, 7.3 and other platforms which allow to 
set ROOT package 

Programming language used: C++: gcc 2.96 or 2.95.2 compiler suite with g++ 
Size of the package: 2.3 MB (0.4 MB compressed distribution archive), without 
ROOT hbraries (additional 120 MB) and without input files with events 
Distribution format: tar gzip file 

Additional disk space required: Depends on the number of events: 35 Mb for 100 000 

events (output of SIMUB generator) 

PACS: 02.70.Tt; 02.70.Uu; 07.05.Tt; 13.25.Hw; 

Keywords: particle physics, decay simulation, Monte Carlo methods, exclusive B- 
meson decay, angular moments method, CP-violation 

Nature of the physical problem: The package BtoVVana has been developed to study 
the performance of the angular moments method of the AP determination from anal- 
ysis of untagged decays — > J/'0(— > //"'"/x~)0(— > K'^K"). By using the package 
BtoVVana it was shown that the method of angular moments gives stable results and 
is found to be an efficient and fiexible tool for measurements with B'^ J/i^P and 
Bl J/^K* decays. 

Method of solution: The method of angular moments allows to construct a sequential 
chain of extraction of the physics parameters from events: AP — >■ initial transversity 
amplitudes weak CP- violating phase strong CP-conserving phases AMg. 
Standard angular moments method was improved by using the time dependent 
weighting function and time-integrated obscrvablcs with variable upper time limit. 
These observables allow to keep full time informative contents of the decay. This 
improved method is used in the package BtoVVana. It was shown that the statistical 
errors of extraction of the observables strongly depends on the choice of the angular 
weighting functions. The best weighting functions were found and implemented in 
the package BtoVVana. The program BtoVVana includes methods of extractions of 
observables as independent modules. In frame of the package BtoVVana user can 
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include new analysis modules. The package BtoVVana and program codes for the 
channels — > J/^(f) and J/'tjjK* in the Monte Carlo generator SIMUB were 

written in parallel and mutually tested with high precision. 

Restrictions on the complexity of the problem: Program processes any number of files 
with any number of events; program tested with about 200 files including 2 x 10^ 
events 

Typical running time: On a PC/Linux with 2.0 GHz processors BtoVVana package 
spends 15 seconds for every 100 000 events. 
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1 Introduction 



The package BtoVVana has been developed to analyze the decays of neutral pseu- 
doscalar B-mesons into two vector mesons decaying into two muons and two pseu- 
doscalar mesons. Two channels of this type are included in the B-physics generator 
SIMUB [1] used here for precise testing of the package BtoVVana: _B°(t),5^(t) —>■ 
J/i^{-^ K+K-) and B^^{t) ,'^^{t) J/il^{-^ n+n-)K*{^ Kn). These 

channels contain rich physics information because of nontrivial angular distributions 
and time dependence of the decays. 

At present we have the information from DO experiment at the TEVATRON 
about observation of 337 events with decay of the first type S° — > J/t/jcf) and 1370 
events of the second type J/t/jK* [2]. The DO data sample corresponds to 

an integrated luminosity 0.22 fb~^, collected in 2002-2004. The statistics will be 
increased up to 4000 J/ip(p events at the DO in the nearest future. Assuming 

the bb cross section to be 100 //b, the number of expected B^ — > J/t/jcf) events, using 
2 fb"-^ of integrated lumunosity, will be equal to 41400 in the BTEV experiment [3] . 
About 8500 B^ — > J/ipK* events are reconstructed in a date sample taken by the 
BELLE detector in the KEKB e+e" collider [4]. 

The systematic study of these decays will be performed after 2007 year on the 
LHC detectors. About 200 000 B° J /^xp events at the CMS, 100 000 events at 
the ATLAS and 75 000 events at the LHCb are expected to be obtained during first 
year of low luminosity operation [5] . 

Possibilities of the current experiments for these important decays can be im- 
proved by using optimal methods of the analysis. For large statistics the best method 
is likelihood fit while for small statistics this method becomes unstable for big number 
of unknown physics parameters, and the angular moments method is more preferable 
[6]. The latest method allows one to build a scquencial chain of extraction of the 
physics parameters [1]. The goal of the paper is to formulate an optimal scenario of 
this analysis which can be realized for real data by using the package BtoVVana pre- 
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sented here. Theoretical ideas reahzed in BtoVVana are based on papers [6] and [1]. 



2 Structure of the package 

Files of the package BtoVVana are kept under the BtoVVana directory which contains 
subdirectories mak, src, doc, res: directory src contains the source codes of the 
program; mak contains Makefile, command files for compilation make_release and 
execution run; bin contains the results of the compilation (object files) and the 
executable file (this directory is created authomatically by command make_release) ; 
doc includes documentations; res is the user directory for data and results. 
The structure of the directory tree one can be found in Makefile: 

source_dirs =../mak 
source_dirs+=T_run_Read 
source_dirs+=T_run_Read/T_Accum_Measure 
source_dirs+=T_run_Read/T_Accum_Measure/T_AngMomMethod 
source_dirs+=T_run_Read/T_Accum_Measure/T_AiigMomMethod/T_WeightFunc 
source_dirs+=T_ruii_Read/T_Accum_Measure/T_AngMomMethod/T_StaiidardModel_DG 
source_dirs+=T_run_Read/T_Accum_Measure/T_Ang]yiomMethod/T_WeightFunc\ 
/T.TimelntObs 

source_dirs+=T_run_Read/T_Accuiii_Measure/T_AngMomMethod/T_WeightFunc 

/T_TimeIiitObs/T_Physics 
source_dirs+=T_Utility 

The directory tree reflects a logical structure and dependences of the classes of 
the package. Command make_release processes all * . C files in directories included 
in the list above. As a rule the name of the directory coincides with the name of 
the class which is placed in the directory. The structure of the program is shown in 
Fig.l. 

Class T_run_Read is intended for reading the events from external files. This 
class contains the main loop over events. For mode Mode_Loop = 2 the loop calls 



the virtual dummy method which will be overlapped in the derived class. Class 
T_run_Read contains the methods for regime Mode_Loop = 3 in which the angle and 
time distributions are created (see Fig. 3). 

The daughter class T_Accum_Measure overlaps some dummy methods of the mother 
class T_run_Read and is intended for regime Mode_Loop = 2. The main data mem- 
bers of the class are the so called "containers" which acummulate the information 
event by event, average the observables, and write the results to the final hsting. 

The class T_AngMomMethod carries out the extraction of the observables by us- 
ing the angular moments method. Two "containers" are impemented in the class 
T_Acciim_Measure: both of the "containers" are described as objects of the same class 
T_AngMomMethod with different angular weighting functions. 

The main data members of the class T_AngMoiiiMethod are four objects of the class 
T_TimeIntObs and three objects of the class T_StandardModel_DG. 

The class T_TimeIntObs contains descriptions of the time-integrated observables 
which are obtained by accumulation of the event weights. 

The class T_StandardModel_DG performs calculations of the width F and width 
difference AF by using the found observables in case of the Standard Model prediction 
of the small CP-violating weak phase. 

The weights of the events are found by using the class T_WeightFunc which is 
the mother class for class T_TiineIntObs. 

Extraction of the final observables based on some predefined physics parameters 
with its errors which described in class T_Physics. The object of this class should 
be created in the main program by the user in case of processing the real data. In 
case of processing the data from generator SIMUB one needs to set logical parameter 
bRealData = false in the main function. In this case the program will read the 
generator physics parameters from the first event. 

Directory T_Utility includes some auxiliary functions and classes. 

The structure of the package BtoVVana is the same as the structure of the gener- 
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ator SIMUB [1] which allows to obtain the events B'^{t),B^{t) J/ip{-^ l+l-)(f){^ 
K+K-) and B^(t),^a(t) ^ J/i^i-^ I'^l') Kn) with full physics contents and 

high precision of generation of the kinematical variables. 

3 Preparation of input files by the SIMUB gen- 
erator 

Subpackage SIMUB/BB_dec allows to obtain the events with Eg — > J/V'^ and — > 
J/t/jK* decays with full physics contents by setting the channel option "COPT SIl" 
in command file BB_dec/mak/run [1] and using the decay mechanism option "BOSM 
1" ("BODM 1") and decay channel option "BOSC 1" ("BODC 1") for 5° (B^) mesons. 
The format of the SIMUB event is defined by option FRMT. Subpackage SIMUB/BB_dec 
has a mode Mode_Loop==2 to produce events in the format readable by the program 
BtoVVana. The variable Mode_Loop is set in the file BB_dec/mak/run by option MODD 
2. In this case the format defined by option FRMT is ignored and the events are written 
in ROOT tree in the format defined in the constructor of class T_Loop (subpackage 
SIMUB/BB_dec) in the following way: 

if (Mode_Loop==2){ 

hfile= new TFile (f name, "RECREATE" , "BOs->J/PsiPhi or BOd->J/PsiK* tree"); 

tree = new TTree("T" , "BOs->J/PsiPhi decay tree"); 

tree->Branch("VB" ,&VB , "x/D:y :z :t :tau") ; 

tree->Branch("PB" ,&PB , "x/D:y:z:E:M") ; 

tree->Branch("Pa" ,&Pa , "x/D:y:z:E:M") ; 

tree->Branch ( " Pal " , &Pal , "x/D : y : z : E : M" ) ; 

tr ee->Br anch ( " Pa2 " , &Pa2 , " x/D : y : z : E : M " ) ; 

tree->Branch("Pb" ,&Pb , "x/D:y:z:E:M") ; 

tree->Branch ( "Pbl " , &Pbl , "x/D : y : z : E : M" ) ; 

tree->Branch ( "Pb2 " , &Pb2 , "x/D : y : z : E : M" ) ; 

tree->Branch( "DFG" , &DFG , "deltal/D : delta2 : Wphi : nA02 : nAP2 : nAT2 : 
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GB:GL:GH:DG:MB:ML:MH:DM:dx:dy:dz:dt:Tmax:cThbl:cThal:sChi: 
cChi : Chi : Blif etime : f Dist : KFCode/I" ) ; } 

Here the f name is a name of the output file defined by option OUTF ' NAME . root ' , 
the branches VB and PB define the secondary 4- vertex of S-decay and 4-momentum 
of B-meson, branches Pa and Pb define 4-momentum of J/ip and {K*) mesons, 
branches Pal and Pa2 define 4-momcntum of ji^ and /i~ mesons, branches Pbl and 
Pb2 define 4-momentum of and K~ mesons {K and tt for K*), the branch DFG 
defines the physics parameters of the decay S — > J/'4>(t> (-6° — > J/'4>K*): strong CP- 
conserving phases ^1 = arg(AjjA_|_) (named as deltal) and 82 = aTg{A±) (delta2; 
we set arg(74o) = 0), weak CP- violating phase (Wphi), initial transversity am- 
plitudes squared A^, A'i, Ay (names nA02, nAP2, nAT2), S° (B^) width F, widths 
of hght and heavy states F^, Fjj, AF = Th — ^l (with the corresponding program 
names GB, GL, GH, DG), B mass M^o and masses of light and heavy states M^, Mh, 
AM = — Mi (program names MB , ML, MH, DM) , generator resolutions [1] or de- 
tector resolution named by dx , dy , dz , dt for four kinematical variables cos 0^+ , 
cos0/^+, X and t, correspondently, where three angles are defined in Fig. 2, and t is 
the B meson proper lifetime. 

If the program BtoVVana is used for the generator data, then the data members 
from the branch DFG contain physics parameters from the generator. In case of 
analysis of real data these data members should contain theoretical predictions to 
calculate some values which depend on the measured values and external parameters 
and for comparison. 

Any number of the input files with events can be treated by the program BtoVVcina. 
The names of these files are written in the main program. 
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4 Quick start with the package BtoVVana 

The package is tested on Linux (RedHat 7.x, 8.x) platforms and uses the ROOT 
package [7]. To run the package, it is necessary to set the environment variables of 
the ROOT. 

To install the program, it is enough to untar the file BtoVVana. tar .gz in the 
working directory and the BtoVVana directory with the structure described above 
will be created. 

The main function is placed in the file inak/T_Accum_Measure_main . C. The user 
can set here the following parameters: 

• code of B-meson KF_Code (Pythia KF-code is 531 for B° and 511 for B^); 

• number of the input files NFiles and its names collected in the array cFile; 

• Mode_Loop = 2 to extract the physics parameters or Mode_Loop = 3 to draw 
the angular and proper time distributions; 

• maximal proper lifetime used to select the events (or for Monte Carlo generation 
in generator) Time_maximal_imnDc; 

• the upper limit for the proper lifetime of B-meson Time_UP_mmDc (in units of 
[mm/c]), Time_UP_mmDc = Time_maximal_mmDc as a rule; 

• second choice of the upper limits for proper lifetime Time_Zero_mmDc; 

• parameter GammaSjnmDcml may be set as arbitrary; to minimize the errors of 
measurements, one needs to set GammaS_mmDcml closely to the true value of F 
(F is B-meson width in units [mm/c]~^); 

• logical parameter bRealData = true in case of analysis of real data, otherwise 
the program will process the events from the generator; the only difference 
is that in case of bRealData = false the user does not need to define the 
object of T_Physics *Phys in the main finction because this object will be 
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defined within the program as the data member by method Paraiiieter_SIMUB 
of the class T_AngMoinMethod by using the special physics records in the events 
written by generator SIMUB (see branch DFG in the class T_run_Read which one 
can found in package BtoVVana and in generator SIMUB); in case of bRealData 
= true the user needs to set an object of the class T_Physics in the main 
function by using the data from the previous run of the program or from other 
sources; 

• number of events NEVENTS (will be ignored if it is less than the number of the 
events in all the input files); 

• to draw a simple histogram, one has to set bSIMPLE_HISTOGR = true otherwise 
it will be drawn histograms with error bars as it is shown in Fig. 3; 

• logical variable DEBUG = true for extended information in the listing used for 
debug of the program; 

• logical variable bShortINF = true for short information in the listing; 

To compile the program, one has to go into mak directory and execute the com- 
mand maJ5e_release. After that, one should start the program by the command 
run. 

Example of the main program and listing with results of the program are placed 
in section " Test Run Input and Output" . 

5 Class T_run_Read: loop over the events 

The main part of the class T_run_Read was done authomatically by the ROOT method 
MakeClass of the TTree class. The constructor of the class T_run_Read contains the 
following parameters: KF_Bmeson is a PYTHIA KF-code of B-meson, NumbFiles is 
a number of the files for analysis, Mode_Loop is a processing mode, and cFile [] is a 
hst of the file names. 
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Class T_run_Read contains the loop which reads the data and writes the event 
information in the folowing data members of the class: 

double VB_x, VB_y, VB_z, VB_t, VB_tau, // Vertex and lifetime 



PB. 


-X, 


PB. 


-y. 


PB. 


-Z, 


PB_E, 


PB_M, 


// 


B momentum and mass 


Pa. 


-X, 


Pa. 


-y. 


Pa. 


-Z, 


Pa_E, 


Pa_M, 


//J/psi moment™ and mass 


Pal. 


-X, 


Pal. 


-y. 


Pal. 


-Z, 


Pal_E, 


Pal_M, 


// 


mu+ momentum and mass 


Pa2. 


-X, 


Pa2. 


-y. 


Pa2. 


-Z, 


Pa2_E, 


Pa2_M, 


// 


mu- momentum and mass 


Pb. 


-X, 


Pb. 


-y. 


Pb. 


-Z, 


Pb_E, 


Pb_M, 


// 


Phi/Ktar momentum and mass 


Pbl. 


-X, 


Pbl. 


-y. 


Pbl. 


-Z, 


Pbl_E, 


Pbl_M, 


// 


K+(K) momentum and mass 


Pb2. 


-X, 


Pb2. 


-y. 


Pb2. 


-Z, 


Pb2_E, 


Pb2_M; 


// 


K-(pi) momentum and mass 



// Physics, kinematics and Monte Carlo parameters (branch DFG) : 
double DFG_deltal ,DFG_delta2,// strong CP conserving phases 
DFG_Wphi, // weak CP violating phase 

DFG_nA02,DFG_nAP2,DFG_nAT2, // initial polarized amplitudes squared (0, I I ,T) 
DFG_GB ,DFG_GL,DFG_GH,DFG_DG,// see definition in previous section 
DFG_MB,DFG_ML,DFG_MH,DFG_DM, // see definition in previous section 
DFG_dx,DFG_dy,DFG_dz,DFG_dt, // see definition in previous section 
DFG_Tmax, // maximal alowed B lifetime 

DFG_cThbl,DFG_cThal, // cos (Theta_K+) cos (Theta_mu+) 

DFG_sChi, DFG_cChi,DFG_Chi, // cos (Chi) , sin(Chi) , Chi 
DFG_Blif etime, // B Lifetime 

DFG_fDist; // value of distribution function for given 

// cos (Theta_K+) , cos (Theta_mu+) , Chi, B-lif etime 
int DFG_KFCode; // PYTHIA KF code of B meson 

Using the method Moiii_To_LorentzVector the event momenta and vertex are 
rewritten in ROOT Lorentz vector variables: P_B_L, P_a_L,P_al_L, P_a2_L, P_b_L, 
P_bl_L, P_b2_L, V_B_gen, V_B_dec, DVertex_B = V_B_dec - V_B_gen. The data 
members - cThe_bl_R, cThe_al_R, Chi_R, B_lif etime_R - present the values cos 6^^+, 
cos©x+, X, t (see Fig. 2) which are calculated for the current event by the method 
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Reco_Helicity_Angles (if the position of a method is not given, then the method 
belongs to the class T_run_Read in this section by default). 

The number of the input files for the processing is set in the constructor of the 
class T_run_Read. by parameter N_umbFiles. The file names are collected in the 
character array in the main program and passed in the constructor by the parameter 
cFile[]. 

The method Loop(NEVENTS) performs the loop over all the events in the file 
chain. The user can find two logical keys bSIMUB_Angles and bSIMPLE_HISTOGR 
inside the method Loop. If bSIMUB_Angles = true, then the angles and B hfetime 
directly from generator SIMUB (from the branch DFG) will be used for histograms, 
otherwise the angles and B lifetime which were reconstructed from the momenta and 
vertex, are used for histograms; if bSIMPLE_HISTOGR = true, then we can obtain 
the pictures with simple histograms for the three angles and B lifetime, while for 
bSIMPLE_HISTOGR = false one can obtain histograms with error bars. The latest 
kind of pictures for decay —>■ J/ip4> can be seen in Fig. 3 for 100 000 events. We 
have used the "main setting" of physics parameters by default for S° — > J/t/j(f){B^ — > 
J/i/jK*) decays: 5i = tt, ^2 = 0, (/)(^) = 0.04, Al = 0.54(0.56), Al = 0.16, Af^ = 
l-Al-Al,T = I/tb, tb = 1.464(1.571) ps, AT/T = -0.2(-0.01), T„,, = 2mm/c, 

Nreso = 50 000, dt^dx^dy^ 'l/Nreso. dz = 27r /Nreso, T^aa^ ^ 2 ^ , X = AM/F = 

20(0.73), where the values for — > J/t/jK* are shown in brackets if they do not 
coincide. 

Switch of regimes Mode_Loop are placed in the main program. Class T_run_Read 
includes methods in a simple mode Mode_Loop = 3. Methods for the Mode_Loop = 
2 are placed in the derived class T_Accum_Measure (see the next section). All other 
values of Mode_Loop arc not used in the current version of the program. 

It is convenient to create a separate derived class TJVccumJIeasure specially to 
treat the events. 

The main loop over the events is placed in class T_run_Read with non virtual 
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definition and therefore it is not redefined in tlie derived class T_Accum_Measure. 
The program uses the same method Loop placed in the parent class T_run_Read for 
all regimes. 

The method Loop contains the switch on different regimes of processing by the 
key Mode_Loop. This switch has the following form (it is a simplifyed form of the 
loop) : 

for(jeiitry=0; j entry<NEVENTS ; jentry++) { 

if (Mode_Loop==2) 
if ( ! Reco_Acciain_Measure ) out("$WARNING Kgoto next evt)",0); 

else if (Mode_Loop==3) Fill_Hist_RECO() ; } 
if (Mode_Loop==2) Average_Measure_OutResult () ; 

else if (Mode_Loop==3) Draw_Hist (0) ; 

Within the class T_run_Read there are dummy virtual methods Reco_Accum_Measure 
and Average_Measure_OutResult which are overlapped in the derived class. The 
method Reco_Accuin_Measure transforms the momenta of the particles into angles 
and accumulats observables. The method Average_Measure_OutResult performs 
averaging and printing the results of the observables extraction. 



6 Class T_Accum_Measure: accumulation of the 
observables 

Class T_Accuin_Measure has T_run_Read as a parent class and its constructor has the 
form: 

T_Accum_Measure(bool DEBUG, 

int KF_B_To_Analyse , int NumbFiles, char *cFile [] , 

int Mode_Loop, double TimeUP_iimiDc , double Tiiiie_maxiinal_iiunDc , 

double Time_Zero_mmDc , double GainmaS_mmDcml , 

bool bRealData, T_Physics *Phys) ; 
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where the parameters DEBUG, KFB_To_Analyse, NumbFiles, cFile [] and ModeLoop 
are used for the constructor of the parent class T_run_Read described in the previous 
section. All other parameters are described in section 4. 

At present the class T_Accum_Measure includes the extraction of observables by 
the angular moments method described in [1]. 

The main methods of the class overlap the parent virtual methods Reco_Accum_Measure 
and Average _Measure_OutResult. These methods are called by derived method 
T_run_Read: :Loop which is not overlapped in the class T_Accuin_Measure. Class 
T_Accum_Measure has the "containers" to accumulate the information in the loop 
over the events. These data are processed and printed in their final form in function 
Average_Measure_OutResult. "Containers" are represented as the following data 
members of the class: 

T_AngMomMethod *AngMomMethod_A ; 

T_AngMomMethod *AngMomMethod_B ; 
where AngMomMethod_A(AngMomMethod_B) are used to treat the events by using the 

angular moments method with the weighting function which is described in [1] as 

"set A" ("set B"). 

The user can introduce new "containers" for other ways of treatment. For ex- 
ample, tagged events of decay Bg{t)B^g{t) — > J/ip(f> containing the oscillation phe- 
nomenon allow to extract the parameter AM. According to the scenario presented 
here, the extraction of the AM by using the tagged samples should be done after 
analysis of the untagged sample. The method to extract AM and corresponding 
"container" will be developed in the next version of the program BtoVVana. 

A pointer to the object of the class T_Accuin_Measure *AM is defined in the main 
function. Constructor of the class T_Accum_Measure defines the "containers" (pa- 
rameters of the constructor will be described in the following section) : 

AngMoiiiMethod_A = new T_AngMoiiiMethod(l , ...); 
AngMoiiiMethod_B = new T_Ang]yioiiiMethod(2, ...); 



13 



Then the main program calls the loop method inherited from the parent class 
AM — > Loop (500000). It passes the number of events to the loop to be processed. 
The events from the external files are read and written into the data members of 
the parent class T_run_Read. The structure of the loop was shown in the previous 
section. As you can see, the loop calls the method Reco_Accum_Measure which is 
overlapped in the class T_Accum_Measure and has the following form: 

bool T_Accuin_Measure: :Reco_Accimi_Measure() { 
Mom_To_LorentzVector () ; 
Reco_Helicity_Aiigles() ; 

Count _MeasureEv++; if (Count _MeasureEv==l) Measure_Init() ; 
Measure_Accuin() ; return true; } 

The first method Mom_To_LorentzVector transforms the momenta to the ROOT 
Lorentz vectors and to the angles (in the frame shown in Fig. 2) and the B-meson 
proper lifetime calculated by using the function Reco_Helicity_Angles. The last 
function is placed in directory T_Utility/F_Reco_HelicityAngles. 

Before treatment of the first event the initialization is performed by means of 
the method Measure_Init. In the treatment of events from the generator SIMUB the 
method Measure_Init calls the method Parameter_SIMUB to initialize the "contain- 
ers" AngMoinMethod_A and AngMoiiiMethod_B by the physics parameters (in this section 
the methods without references have the placement in the class T_Accum_Measure by 
default). 

These parameters are used to extract the observables from the time-integrated 
observables which have been obtained by class T_TiineIntObs, and to perform ana- 
lytical calculations to compare them with the approximate estimations obtained by 
the Monte Carlo method (see description of the class T_Physics below). It is helpful 
to test the both the event generator and the analysing package BtoVVana. 

As it is described in the previous section, the method Loop derived from class 
T_run_Read calls the virtual method Reco_Accuin_Measure which overlapped in the 
class described in this section. Further the chain of calls is the following: Reco_Accuin_Measure () 
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Measure_Accuin() — > Measure_Accuin_AngMomMethod() . The structure of the lat- 
est method is very generah 

void T_Accum_Measure : :Measure_Acc;im_AngMomMethod(){ 
AngMoinMethod_A->Init_Angles(cThe_bl_R, cThe_al_R, Chi_R, B_lif etime_R) ; 
AngMoinMethod_A->Accuiiiulation_TimeIntObs () ; 
AngMoinMethod_A->Set_False_FlagInit_Angles() ; 

AngMomMethod_B->Init_Angles(cThe_bl_R, cThe_al_R, Chi_R, B_lif etime_R) ; 
AngMoinMethod_B->Accuiiiulation_TimeIntObs () ; 
AngMoinMethod_B->Set_False_FlagInit_Angles ( ) ; 

} 

It is not difficult to include other user "containers" into this method. 

In method T_run_Read: :Loop(int) after the loop the method 
Average_Measure_OutResult is called. This method performs averaging and print- 
ing the results of the measurements (it is a simplified form) : 

void T_Accum_Measure : :Average_Measure_OutResult(){ 
AngMomMethod_A->Average_TimeIntObs () ; 
AiigMomMethod_A->Out_Result_TimeIiitObs(" Set A "); 
AngMoinMethod_B->Average_TimeIntObs () ; 
AngMoinMethod_B->Out_Result_TimeIntObs(" Set B ");> 



7 Class T_AngMomMethod: extraction of physics 
parameters by using the angular moments method 

The main data members of the class T_AngMomMethod are four pointers to the objects 

of the class T_TimeIntObs and three pointers to the objects of the class T_StandardModel_DG: 



T. 


.TimelntObs 


*Obs. 


.TGO ; 


// 


T , 


Gp = 





T. 


.TimelntObs 


*Obs. 


-TGp ; 


// 


T , 


Gp = 


G' 


T. 


.TimelntObs 


*Obs. 


.TOGO; 


// 


TO, 


Gp = 
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T.TimelntObs *Obs_TOGp; // TO, Gp = G' 
T_StandardModel_DG *SM_DG_TGp_TGO ; 
T_StandardModel_DG *SM_DG_TGp_TOGp ; 
T_StandardModel_DG *SM_DG_TGO_TOGO ; 

Two different definitions of the time-integrated observables can be considered 
according to [1] : O and O. The first case corresponds to the zero parameters F' 
and the second one - to nonzero V defined by the user (see the next section about 
definition of F'). In the program the user can set also two different values of the 
upper time limit T and To- 

Constructor and methods of the class T_AngMoiiiMethod are 

T_AngMomMethod (...); 

void Init_Angles (...); 

void Parameter.SIMUB ( . . .) 

void Accumulatioii_TimeIntObs () ; 

void Average_TimeIntObs () ; 

void Out_Result_TimeIntObs (const Char_t *cc) ; 

virtual void Show (const Char_t *c) ; 

The constructor calls four constructors for four objects of the class T_TimeIntObs. 
Each of the methods listed above calls four methods with the same name for the four 
objects of the class T_TimeIntObs and for the three objects of the class T_StandardModel_DG. 
For example, 

void T_AngMomMethod : :Acciunulation_TimeIntObs(){ 
Obs_TGO ->Acciamulation_TimeIntObs() ; 
Obs_TGp ->Accuinulation_TimeIntObs ; 
Obs_TOGO->Acciainulation_TiineIntObs () ; 
Obs_TOGp->Accuiiiulation_TimeIntObs () ; 
SM_DG_TGp_TGO ->Accumulation() ; 
SM_DG_TGp_TOGp->Accumulation() ; 
SM_DG_TGO_TOGO->Accumulation() ; > 
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Some features in using the objects of the classes T_TimeIntObs and T_StandardModel_DG 
the user can find in the method Parameter_SIMUB (this method calls the methods of 
T_TimeIntObs objects only), and in the method Average_TimeIntObs: 

void T_AiigMomMethod : :Average_TimeIiitObs(){ 
Obs_TGO ->Average_TimeIntObs() ; . . . ; 

SM_DG_TGp_TGO->T_StaiidardModel_DG_Init(*Obs_TGp, *Obs_TGO ); 
SM_DG_TGp_TGO->Average() ; SM_DG_TGp_TGO->Calc_DG(true) ; . . . ;} 

Method T_StandardModel_DG_Init initiates the object to extract F and AF. 
Method Average averages values to calculate the correlation between two types of 
the observables. Method Calc_DG calculates F and AF. 

8 Class T_WeightFunc: weighting functions 

Two types of the angular weighting function are proposed in [6, 1] to extract the 
observables: wI'^\6i+,6k+,X) ^^'^ w\^\0i+^9k+-,x) = 1)6). As it is shown in [1], 
it is helpful to consider the time dependent weighting functions which can be written 
in the general case as 

wi^\ei+, 9k+,x, t; r', T) = eM^'t)w\^\ei+,eK+.xmT - t) 

and similarly for H^/'^^ . Weighting function wj;^^^'^ is named in the class T_WeightFunc 
as method wi {i = 0, 5). 

The integer data member WF_Type of the class T_WeightFunc sets the type of 
weighting functions {1 —>■ A, 2 ^ B) and it is initialized by the class constructor. The 
parameters of the weighting functions F' and T are the constructor input parameters, 
named as GanraiaS_mmDcml (in unit (mm/c)~^) and TimeUP_minDc (in unit (mm/c)). 
The F' represents a first approximation for the measured value F to be corrected in 
the current analysis of the untagged decays. 

Together with method wi there are also the methods with names wi_cThal, 
wi_cThbl, wi_Chi and wi_t in the class T_WeightFunc. These methods represent 
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the derivatives of the function wi with respect to arguments cos0^+ , cos0x+ , X and 
t, correspondently. The derivatives are used in class T_TimeIntObs to calculate the 
systematical errors or to estimate the detector response [1] . 



9 Class T_TimeIntObs: six time-integrated ob- 
servables 

Class T_TiineIntObs is a derived class to the class T_WeightFunc. It is used to 
accumulate and average the time-integrated observable, and to extract the physics 
parameters. 

In the general form by means of the weighting function wj;^'^ (X — A, B) on 
the set of the N{Tmax) events we have observables and their statistical errors in the 
following form [1]: 



r^^(r', T) = — ^ E wP{T\ n (1) 

' I,-' max ) j=i 



N{T) 



N{T) 



i=i 

where W^piV, T) = W^^^\ej+,e^j^+,x^ , V; T', T) is a value of the weighting function 
for the j-th event. Including the time dependent 0(T — i)-function in the weighting 
function allows to keep the time informative contents in the value })f '^^^(T\ T) and 
in the same time to have good statistical errors which are caused by a large size of 
the proper time interval [0, T\. 

The class constructor T_TimeIntObs is defined by the following signature: 

T_TimeIntObs(char *cc, int Type_WF, 

double GaiimaS_iiiiiiDcml , TiiiieUP_iiiiiiDc , TimeMAX.mmDc , 
bool DEBUG, bRealData, T_Physics *p_hys) ; 

Three parameters Type_WF , GaimnaSjimDcml , TimeUPjnmDc are sent to the mother 
class T_WeightFunc (see the previous section). Parameter TimeMAXjmnDc defines the 
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By means of the logical parameter bRealData the user can set either the mode 
for the real data or the mode for the generator data. The latest case does not use 
the method parameter P_hys and data member Phys of the class will be filled by the 
generator data taken from the first event by the method Parameter_SIMUB. 

In case of bRealData = true the user has to define the object of the class 
T_Physics in the main program and send it into the object of the class T_Accuin_Measure 
defined in the main program. 

The constructor and hst of the methods of the class are as follows: 

T_TimeIiitObs(char *cc, int Type_WF, 

double GaiimaS_iiiiiiDcml , TimeUP_iiiiiiDc , TimeMAX_iiiiiiDc , 
bool DEBUG, b_RealData, T_Physics *P_hys) ; 
void FinalObserv_TimeIntObs (const Char_t *c) ; 
void Parameter_SIMUB (double deltal, delta2, TotalPhaseWeak, 

nA0_02, nA0_P2, nA0_T2, Gamma_Biiies_GeV, Gp_GeV, Giii_GeV, Delta_G_GeV, 

Gen_dcTbl, Gen_dcTal, Gen_dChi, Gen_dt) ; 
void Accumulation_TimeIntObs() ; 
void Average_TimeIntObs () ; 
void Out_Result_TinieIntObs (const Char_t *c) ; 
virtual void Show (const Char_t *c) ; 

If we compare this list with the list of the methods of the class T_AngMomMethod, 
it is possible to note that the last 5 methods have the same names as the method 
names of the class T_AngMomMethod because the calls of these methods are the main 
purpose of the corresponding methods of the class T_AngMoinMethod (see explanation 
in the description of the class T_AngMomMethod). 

Method FinalObserv_TimeIntObs is called by method Out_Result_TimeIntObs 
and calculates the final observables or their combinations according to the formu- 
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lae [1]: 



Ur''\v\ T) + &f"^)(r', T) + ht''\V\ T)/7 



|^,(0)P- r'>(r,r)/, 



T) + 6f"^)(r', T) + 6f"^)(r', T)/7 ' 

oos{5, - 5,) = ^ : \ , (2) 

where we consider the initial amphtudes normahzed as | Ao(0) P + | Ay (0) p + |^±(0) p 
= 1. We have introduced the function 7: 

7(Ari,Ar^,cos0W,T) = ^, (3) 

G . = (1 ± cos0(^)) (IT cos</.i^)) . (4) 

AFl/h are measured parameters related with physics values of widths of light and 

heavy states Ti/h by 

Ar^ = 2(r' - r^) , Ar^^ = -2(r' - r^^) . (5) 

The method FinalObserv_TimeIntObs also calculates two values: 

sin0(^) cos5,,. = '''' p. (6) 

where 



/?(A^i,A^^^,cos0^^T) = ^^^^ , Z= ^^^^ + . (7) 

Eqs. (6) contain the weak phase in the left- and right-handed side and might be 
helpful in case of large values of weak phase (p^J'^ when a large violation of the Standard 
Model predictions, 0^*) ^ 0.03, takes place. 

Values 7 and $ are calculated in the data member T_Physics *Phys of the class 
T_TimeIntObs. 
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Table 1: Time-integrated observables obtained by data member Obs_TGO of class 
T_AngMoinMethod. The sample of 300 000 (3 000 in parentheses) untagged S° J/ipcf) 
events from SIMUB generator with the "main setting" of physics parameters (see 
section 5) is used 



i 








Statistical 




Eq. (8) 


Eq. (8) 


Eq. (1) 


errors 


1 


0.54 


0.5225 


0.5207 (0.488) 


0.0014 (0.014) 


2 


0.30 


0.2903 


0.2940 (0.380) 


0.0021 (0.021) 


3 


0.16 


0.1873 


0.1849 (0.132) 


0.0020 (0.020) 


4 


0.0088 


0.0009 


0.0032 (0.011) 


0.0020 (0.020) 


5 


-0.4025 


-0.3894 


-0.3916 (-0.316) 


0.0029 (0.029) 


6 


-0.0118 


-0.0012 


-0.0021 (-0.002) 


0.0031 (0.031) 



The results of extraction of final observables (2) and (6) are given in Tabls. 1 
and 2. 

The direct numerical calculations have shown that the difference between the 
values of observables hi{T) {i = 1,2,3,5), calculated with cjy'f^ = and = 0.04, 
does not exceed 0.01%. Even in case of statistics of 100 000 events this difference is 
negligibly small as compared with statistical errors for these observables. Therefore, 
the assumption 0^*^ f^i is a good approximation in case of Standard Model. 

To determine the initial transversity amplitudes in case of SM, it is convenient 
to use r' = and for large T we have 7 fti Tl/Th- 

The first observable in the fist of four observables of the class T_AngMoinMethod 
is Obs_TGO which is initialized with setting F' = and with T = T^ax as the upper 
time limit. The results from Obs_TGO arc presented by two tables in the listing (see 
section "Test Run Input and Output"). 

These data are given in Tabls. 1 and 2. 
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Tabic 2: Determination of values (2) and (6) by using the observables bf'^-^\o,T) 
{T = 2/1/0.5/0.25/0.125 mm/c"^) for 300000 5° J/il;(f) events with the "main 
setting" (see section 5) of physics parameters 



Parameter 


Input 


Extracted value 


Statistical 






by (2), (6) 


error 




0.556 


0.565/0.568/0.565/0.568 


0.005/0.005/0.006/0.007 


AllAl 


0.296 


0.294/0.293/0.293/0.288 


0.004/0.004/0.005/0.006 


IA(o)P 


0.54 


0.538/0.538/0.538/0.539 


0.001/0.001/0.002/0.002 




0.30 


0.304/0.305/0.304/0.306 


0.002/0.002/0.003/0.003 




0.16 


0.158/0.157/0.158/0.155 


0.002/0.002/0.002/0.003 


cos(52 - 5i) 


-1 


-l,001/-l,000/-0.998/-0.999 


0.009/0.009/0.010/0.012 


cos(5isin0^*') 


0.04 


0.15/0. 17/0.16/-0.04 


0.09/0.12/0.22/0.50 


cos52sin0^'*^ 


-0.04 


-0.07/-0.07/0.03/-0.03 


0.10 /0.14/0.26/0.61 
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Table 3: Time-integrated observables obtained by data member Obs_TGp of class 
T_AngMoinMethod for 300 000 (3 000 in parentheses) untagged S° J/il^(t> events 
with the "main setting" (see section 5) of physics parameters 



i 


6,(r,T^,,), Eq. (8) 


y~r\v, T^a.), Eq. (1) 


Stat. err. 


1 


2.117 


2.139 (1,95) 


0.014 (0.12) 


2 


1.176 


1.157 (1,40) 


0.019 (0.17) 


3 


0.989 


1.024 (0.78) 


0.020 (0.16) 


4 


0.010 


0.035 (0.07) 


0.018 (0.19) 


5 


-1.578 


-1.603 (-0.94) 


0.028 (0.24) 


6 


-0.013 


-0.010 (0.11) 


0.028 (0.31) 



As it is shown in the first two columns of Tabl. 1, the time dependence of ob- 
servables is essential in case the number of events is more than 3 000 because for 
the 3 000 events the differences bi{0,T)T=Tmax^o — bi{0,T)T=Tmax=2 are comparable 
with statistical errors. This conclusion depends on the value of width difference 
Ar. In case of the — > J/t/jK* channel the Standard Model predicts a small 
value of Ar ("main setting" for 5° AT /T — —0.01, see section 5) and in this case 

bi{0,T)T=Tma.-^0 ~ &i(0,T)T=T„„,=2. 

The second observable in class T_AngMoinMethod is Obs_TGp which is considered 
here with parameter F' = F (F is a true value of the B meson width) and with Tmax 
as the upper lifetime limit. The results from listing tables of observable Obs_TGp are 
shown here in Tabl. 3 and 4 and in listing of section "Test Run Input and Output" . 

It is better to use observables b\^^^\o,Tmax) to extract the initial observables (2) 
and (6) (compare the statistical errors in Tabl. 2 and 4). 

Extraction of combinations cos5i^2sin0^*'' is shown in Tabl. 5. 
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Table 4: Determination of values (2) and (6) by using the observables bf^^\r, T^ax) 
for 300000 S° — > events with the "main setting" (see section 5) of physics 

parameters 



Parameter 


Input value 


Extracted vahic 


Stat. error 




0.556 


0.5409 


0.010 


Al/Al 


0.296 


0.3038 


0.0066 


l^(0)P 


0.54 


0.5421 


0.0030 




0.30 


0.2932 


0.0045 




0.16 


0.1647 


0.0034 


C0S((52 - 5i) 


-1 


-1,019 


0.020 


cos(5isin0^''^ 


0.04 


0.138 


0.073 


cos52sin0^*) 


-0.04 


-0.029 


0.081 



Table 5: Determination of values depended on weak phase by using the observables 
^(cxp)^p^ T^„3.) extracted from Monte Carlo data for two values of 0^*-* = 0.4 and 0.04. 
The events sample has been generated for the case of "main setting" for physics 
parameters 



Parameter 


Input value 


100 000 events 


200000 events 


400000 events 


cos(5isin0(*) 


0.04 


0.24 ±0.14 


0.13 ±0.09 


0.12 ±0.07 


cos52sin0^'') 


-0.04 


-0.05 ±0.15 


0.10 ±0.11 


-0.03 ±0.07 


cos5isin0(*) 


0.389 


0.49 ±0.16 


0.44 ±0.11 


0.35 ± 0.08 


cos52sin0^*) 


-0.389 


-0.29 ±0.16 


-0.31 ±0.12 


-0.41 ±0.09 
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10 Class T_Physics: physics values 

For observables 6i(F', T) we have the following formulae [1]: 



61 (F', T) = 


|Ao(0)PGL(r', T)/LiTmax), 




H^', T) = 


|A||(0)|^G'l(F', T)/L{Tma.), 




^3(r', T) = 


\A40)fGH{T', T)/L{Tma.), 




&4(F', T) = 


\A\\{0) \ \A^{0)\Z{r', T)cos5isin 


.(t>^^/L{T,,a.) , 


UT\ T) = 


|Ao(0)||A||(0)|Gi(r, T)cos(52- 


- Si) / L(Tmax) 


T) = 


|Ao(0)| \A^{0)\Z{r', T) cos(52sin 





(8) 

where L{Tmax) is a normalization factor, which has the form: 

L{Tr^ax) = (1^(0) I' +|^||(0)|2) 6-^(0, r„„,) + 1^^(0)12 6^(0, r„„,), (9) 

From expresions (8) one can see that 61 (F', T) + 62(1", T) + 63(F', T) = 1. 

In the general case with sizable weak phase 0^^^ to extract the physics values we 
can use the following equations: 

6S(r,T) 



6(7) (F", To) 



/.^(AFi, AFh, cos0(^); [F', T], [F", To]); (10) 
/iH(AFi, AFh,cos0«; [F',r], [F",To]) (11) 
p(AF,„AFh;[F',T],[F",To]), (12) 



where 



GL//f(AF^,AF^,cos0W;[r,r]) 
G'i/^,(AFi,AF^,cos#;[r",ro])' 
Z(AF^,AF^^;[F',r]) 



(14) 



Z(AF^,AF^;[F",ro]) 
Therefore, to extract tree parameters AF^, AFj^ and cos0^^^ we need the experimen- 
tal values of the time- integrated observables bj^^^^ only 
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In case of small values b^^Q we have a small sin0^*^ cos5i,2, but it does not mean 
the small weak phase 0^*^ as it is predicted by the Standard Model. The opposite 
case of the sizable observables Si^g^^ will directly show the signal of the beyond 
Standard Model physics. In the last case we can extract AF^, and AF^^ from the 
two equations (12). If we use F' as the true value of the B meson width (found 
from other sources), then we will have a single unknown value AF^ = AFj^ = AF. 
Simphfication of the system of equations (10)- (12) in case of the Standard Model 
will be considered below. 

After extraction of AF^, ATh and cos0^*) we can calculate the values 7 (3) 
and /3 (7). The values 7 and $ are calculated in class T_Physics and used in the 
method FinalObserv_TimeIntObs of the class T_TimeIntObs to calculate the final 
observables (2) and (6). It is the first goal of the class T_Physics. 

The second goal of the class T_Physics is to calculate exact theoretical values of 
the observables (8) to compare them with the approximate values obtained by (1) 
in the class T_TimeIntObs. 

For these purposes the class T_Physics collects all the physics values including 
transversity amplitudes with its errors and three user parameters F', T and T^ax- 

In the main method Init of the class T_Physics one can sec the scenario to 
calculate of the physics values described above. The list of the input parameters 
of the method Init includes 8 physics parameters: deltal and delta2 (strong 
CP-conserving phase ^1,2), TotalPhaseWeak (weak CP- violating phase 0^*-*), nA0_02, 
nA0_P2 , nA0_T2 (initial transversity amplitudes squared |^o(0)P, |^||(0)P and |^±(0)p = 
1 — 1^40(0) p — |74||(0)p), Gp and Gm (widths of light and heavy B meson states F^ 
and Th), their errors, and 3 user defined parameters: TimeUP (upper time limit 
T), TimeMAX (maximal time Tmax) and Gprime (F') to use them for calculation by 
formulae (8). Data member oEkJiat {k — 0,5) of the class T_Physics are the 
observables (8) in case of F' defined by user while oEk_tilde {k — 0, 5) are the 
observables (8) in case of F' = 0. 
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11 Class T_StandardModel_DG : extraction of T 
and Ar in case of the Standard Model 



The name of the class T_StandardModel_DG means that it includes the extraction of 
r and Ar in case of the Standard Model expectation for the weak phase: (f)'f^ ^ 0. 
In this case we may not use the small values of the observables 64^5^^ and, according 
to (4) and (7), we have: fi^ = /XilAr^; [r',T], [r',To]), = /ii.(Ar^; [r',T], [r",To]) 
It means that extraction of Ar^ and ATh can be performed by numerical solving of 
the separate equations. 

Three cases numbered as KK = 1,2,3 (see method mu_DG_Equa_Solut) are real- 
ized in the class T_StandardModel_DG: KK = 1 in case of V ^ 0, T" = 0, Tq = T 
(used in data member SM_DG_TGp_TGO of the class T_AngMomMethod) , KK = 2 in case 
of r' = r" = 0, To 7^ 7" (used in data member SM_DG_TGO_TOGO) , and KK = 3 in case 
of r' = r" ^0,To^T (used in SM_DG_TGp_TOGp) . Only the case KK = 1 is presented 
in listing which is shown in the section " Test Run Input and Output" . In each case 
we have 3 equations (10) to determine AF^ (names DGL_0, DGL_1 and DGL_4 in output 
listing) and one equation (11) for ATh (with name DGH_2 in output listing). 

Combining three values AF^ from (10) with one value AFj^ from (11), we have 
three value AF = 0.5(AFz, + AF^j) (names DG_0, DG_1 and DG_4 in listing of section 
"Test Run Input and Output" ) and three value F = F' - 0.25(AFl - ATh) (names 
G_0, G_l and G_4 in listing) corresponding to combinations of indices [1,3], [2,3] and 
[5,3] of observables. 

Three cases KK=1,2,3 of extraction of F and AF are presented as three tables in 
listing. 

The results of the first method of extraction of F and AF (KK=l) weakly de- 
pendend on user defined F' value (see Tabl. 6 obtained by SM_DG_TGp_TGO) . Due to 
weak dependences of the result on F' we fix this parameter in the following tables as 
F' = F. 
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Table 6: Determination of values F and Ar by using the object SM_DG_TGp_TGO (KK 
= l) for 300000 — > J/t/jcf) events with the "main setting" of physics parameters. 
Combinations of indices [1,3] [2,3] [5,3] are explained in section 11 



Parameter 


r' 


Input 


[1,3] 


[2,3] 


[5,3] 


r, [mm/c]~^ 

Ar. [mm/c]-i 


2.278 

2.278 


2.278 

-0.156 


2.228 ±0.015 

-0.506 ± 0.030 


2.268 ±0.020 

-0.587 ± 0.010 


2.231 ±0.020 

-0.512 ±0.0 10 


r, [mm/c]""*" 
Ar, [mm/c]^"^ 


2.392 
2.392 


2.278 
-0.456 


2.227 ±0.015 
-0.508 ±0.030 


2.268 ±0.020 
-0.591 ± 0.040 


2.230 ±0.020 
-0.515 ±0.041 


r, [mm/c]~^ 
Ar, [mm/c]-^ 


2.164 
2.164 


2.278 
-0.456 


2.229 ±0.015 
-0.504 ±0.029 


2.268 ±0.019 
-0.582 ± 0.039 


2.232 ± 0.020 
-0.509 ±0.039 



The results of the second method of extraction of r and Ar (KK = 2) depend on 
the user defined time up limits T and Tq (see Tabl. 7 obtained by SM_DG_TGO_TOGO) . 
The best choice is Tq = 0.5 mm/c in case of maximal T = T^ax = 2 mm/c. 

The results of the third method of extraction of r and Ar (KK = 3) depend on 
time upper hmits T and To defined by the user (see Tabl. 8 obtained by SM_DG_TGp_TOGp) . 
The best choice is Tq = 0.25 mm/c in case of maximal T = Tj^^x = 2 mm/c. The 
last method of extraction of r and Ar is better among all the considered methods. 

To avoid the superfluous information, the user needs to set DEBUG = false and 
bShortINF = true. In this case the listing consists of two parts for two methods 
of treatment by using the "set A" and "set B" weighting functions [1] (in section 
" Test Run Input and Output" you can see the hsting with the case of " set B" only) . 
The " set B" weighting functions are the hnear combinations of six angular functions 
which define the amplitude of the process while the " set A" weighting functions are 
not expressed linearly via the angular functions. This is the main reason why the 
statistical errors for observables in case of "set B" is about 2 times smaller than for 
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Table 7: Determination of values T and AT by using the object SM_DG_TGO_TOGO (KK 
= 2) for 300000 — > J/t/jcj) events with the "main setting" of physics parameters. 
Combinations of indices [1,3] [2,3] [5,3] are explained in section 11 



Parameter 


To;T, 


Input 


[1,3] 


[2,3] 


[5,3] 




mm 

c 


[mmj-i 


[mmj-i 


[mmj-i 


[mmj-i 


r 


i;2 


2.278 


2.246 ±0.017 


2.281 ±0.022 


2.255 ± 0.022 


Ar 


i;2 


-0.456 


-0.478 ± 0.033 


-0.549 ±0.043 


-0.498 ± 0.043 


r 


0.5; 2 


2.278 


2.242 ±0.015 


2.272 ±0.019 


2.243 ±0.019 


Ar 


0.5; 2 


-0.456 


-0.485 ± 0.029 


-0.544 ±0.038 


-0.487 ±0.039 


r 


0.25; 2 


2.278 


2.240 ±0.021 


2.249 ± 0.025 


2.242 ± 0.025 


Ar 


0.25; 2 


-0.456 


-0.509 ± 0.042 


-0.527 ±0.049 


-0.512 ±0.050 


r 


0.125; 2 


2.278 


2.257 ±0.028 


2.236 ± 0.032 


2.210 ±0.033 


Ar 


0.125; 2 


-0.456 


-0.477 ±0.056 


-0.434 ±0.065 


-0.383 ± 0.066 
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Table 8: Determination of values T and Ar by using the object SM_DG_TGp_TOGp (KK 
= 3) for 300000 — > J/t/jcj) events with the "main setting" of physics parameters. 
Combinations of indices [1,3] [2,3] [5,3] are explained in section 11 



Parameter 


To;T, 


Input 


[1,3] 


[2,3] 


[5,3] 




mm 

c 


[mmj-i 


[mmj-i 


[mmj-i 


[mmj-i 


r 


i;2 


2.278 


2.225 ±0.019 


2.279 ± 0.026 


2.238 ±0.026 


Ar 


i;2 


-0.456 


-0.514 ±0.039 


-0.622 ±0.053 


-0.541 ± 0.053 


r 


0.5; 2 


2.278 


2.228 ±0.015 


2.268 ± 0.020 


2.231 ± 0.020 


Ar 


0.5; 2 


-0.456 


-0.506 ±0.030 


-0.587 ±0.040 


-0.512 ±0.040 


r 


0.25; 2 


2.278 


2.231 ±0.014 


2.263 ±0.018 


2.237 ±0.018 


Ar 


0.25; 2 


-0.456 


-0.508 ±0.029 


-0.573 ±0.036 


-0.520 ±0.037 


r 


0.125; 2 


2.278 


2.238 ±0.016 


2.256 ±0.019 


2.222 ± 0.020 


Ar 


0.125; 2 


-0.456 


-0.496 ± 0.032 


-0.532 ±0.039 


-0.465 ± 0.039 
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the "set A". In previous sections we have described the results for "set B" only. 

In previous sections we have used three keys for tuning the program BtoVVana: 
r', T and Tq. To improve solving the equations for AF^^i/, the user can change hmits 
of arguments in the method mu_DG_Equa_Solut (class T_StanciardModel_DG) which 
solves the equations. 

12 Conclusion 

Extraction of physics information in decays —>■ J/ip4> and —>■ J/ipK* by using 
the angular moments method with time dependent and time-integrated observables, 
has a number of attractive features, which are demonstrated by package BtoVVana: 

• it is an unbinned method; 

• it uses full time informative contents of time-dependent decays; 

• it uses full available statistics; 

• it gives stable results in case of small statistics; 

• it allows one to separate the extraction of physics values; 

• it allows one to use different scenarios in case of the signals beyond Standard 
Model, or in case when it is justified; 

• it is a flexible tool because it has a different ways to tune the extraction of 
observables (for example, to tune solving the equations) in the process of the 
real data treatment; 

• and it is a visual method (see, for example, the Fig.4 in [1] for dependences of 
61,2 on Ar). 

The detailed tests have been performed for the package BtoVVana by means of a 
precise generator SIMUB. The tests have checked mutually the both programs SIMUB 
and BtoVVana with high precision. 
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The complex BtoVVana - SIMUB can be used to test other methods of extraction 
of physics information from decays — > J/V'^ and — > J/t/jK*. 

The program BtoVVana has clear structure and can be used as template to include 
new methods of treatments. 
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Test Run Input and Output 

User defined parameters are defined in section 4 and placed in main program: 

TROOT root ("Program BtoVVana" , " xxx "); 
int mainCint argc, char **argv) { 

TApplication *theApp= new TApplicationC'App" , feargc, argv) ; 

cout«" ***** mainO Start ***** "«endl; 
//============= user set, Files with events and parameters ===== 

const Int_t NFiles = 10; const Char_t *cFile [NFiles] ; 

cFileE 0] = "$DATl/10000ev_B0sJPsiPhi_l.root"; 

cFileC 1] = "$DATl/10000ev_B0sJPsiPhi_2.root"; 

cFile[ 9] = "$DATl/10000ev_B0sJPsiPhi_10.root"; 

lnt_t KF_Code =531; // 531 -> BOs.BbarOs, 511 -> BOd.BbarOd 

Bool_t DEBUG =false;// switch for extended listing 

Bool_t bShortINF =true; 

Int_t Mode_Loop =2; // 2: DGamma measurement; 3: angle distributions 
Double_t Time_UP_mmDc=2 . , Time_maximal_mmDc=2 . , Time_Zero_mmDc=0 . 5 ; //mm/ c 
Double_t GammaS_nimDcnil ; 

if (KF_Code ==531) GammaS_nmiDcml = 2.278; 

else if (KF_Code==511) GammaS_mmDcml = 2.27844;// Gd = 2.12326; 
else out_exit("main: check KF_Code or set GammaS_mmDcml . STOP."); 
Bool_t bRealData = false; 

Double_t ErrFactor=0 . 2 , Reso_CosThetap=0 . 022 , Reso_CosThetaKp=0 . 0076 , 

Reso_Chi=0.04,Reso_t=0.03; // [Reso_t] = mm/c 
Int_t NEVENTS = 100000; 
Bool_t bSIMPLE_HISTOGR = false; 
//============= end user setting. ================ 

T_Physics *Phys = new T_Physics(); 

Phys->Init(Time_UP_mmDc, Time_maximal_mmDc , 0., ErrFactor*0 . , 
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3.14, ErrFactor*3. 14, 0., ErrFactor*0 . , 0.33, ErrFactor*0.33, 
0.33, ErrFactor*0.33, 0.33, ErrFactor*0.33, 
true , GaimnaS_mmDcml*Phys->mmDcml_To_GeV , 
GainmaS_inmDcml*Phys->inmDcml_To_GeV*1.2, // GL 
ErrFactor*GanmaS_nimDcnil*Phys->inniDciiil_To_GeV*l .2, // GL err 
GaminaS_iimiDcml*Phys->iimiDcml_To_GeV*0.8, // GH 
ErrFactor*GainmaS_mmDcml*Phys->]iiiiiDcml_To_GeV*0.8, // GH err 
Reso_CosThetap, Reso_CosThetaKp, Reso_Chi, Reso_t) ; 
if (Mode_Loop==2) { 
T_Accum_Measure *AM = new T_Accum_Measure (DEBUG, bShortINF, 
KF_Code, NFiles, cFile, Mode.Loop, 
Time_UP_mmDc , Time_maximal_mmDc , Time_Zero_iiiiiiDc , 
GainmaS_mmDcml , bRealData, Phys) ; 
AM->Loop(NEVENTS, bSIMPLE.HISTOGR) ; 

} 

else if (Mode_Loop==3){ 
T_run_Read *SRAM=new T_run_Read(KF_Code, NFiles, cFile,Mode_Loop, DEBUG) ; 
SRAM->Loop(NEVENTS, bSIMPLE_HISTOGR, 0); 

} 

if (Mode_Loop==3) theApp->Run() ; 
cout«" ***** mainO End ***** "«endl; 

} 

The output listing in short format has the form: 

***** mainO Start ***** 

ctor T_Accum_Measure 

weight func [set B] 

ctor T_TimeIntObs: 

Gamma'prime = 4.49511e-13 GeV = 2.278 [inm/c]"-l 

T_imeUP_inmDc = 1.01355e+13 GeV-l = 2 mm/c 

T_imeMAX_iiiinDc= 1.01355e+13 GeV"-l = 2 mm/c 
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[initial condition for other three kind of observables is placed here] 
Inform from T_TimeIntObs : :Parameter_SIMUB in case of bSIMUB_Test: 

Gamma_BOs = 4.49598e-i3 GeV = 2.27844 [mm/c] 

DG =-8.99197e-14 GeV =-0.455689 [mm/c]"-! 

Gp(GL) = 4.945588-13 GeV = 2.50629 [nun/c] 

Gm(GH) = 4.04638e-13 GeV = 2.0506 [mm/c] "-1 

0.5*(Gm+Gp)= 4.49598e-13 GeV = 2.27844 [mm/c]"-! 

(Gm-Gp)= -8.99197e-14 GeV = -0.455689 [mm/c]"-! 

DG/G = -0.2 

RESOLUTION: dcTbl, dcTal, dChi, dt = 4e-05, 4e-05, 0.000125664, 4e-05 



BEGIN: T_Accum_Measure : : Average_Measure_OutResult 
*** [set B] *** T_Accum_Measure :: Average _Measure_OutResult 
NEVENT = 100000 TimeUP = 2 mm/c, GammaS = 2.278 [mm/c] ~-l 

BEGIN TGp -2priine- BEGIN — 



oO= 


2. 


, 1242+/- 


-0. 


,024( 


1, 


.1°/.)+/- 


-0, 


,0000(0%)= 


2. 


,1169 


= oEO 


ol= 


1, 


, 0584+/- 


-0, 


,033( 


3, 


.1%)+/- 


-0, 


,0000(0%)= 


1. 


,1760 


= oEl 


o2= 


1, 


,0787+/- 


-0, 


,035( 


3, 


.2%)+/- 


-0, 


,0000(0%)= 


0. 


,9889 


= oE2 


o3=- 


-0, 


.0165+/- 


-0, 


.031(188.3%)+/- 


-0.0000(0%)= 


0. 


,0099 


= oE3 


o4=- 


-1, 


,5242+/- 


-0, 


,049( 


3, 


,2%)+/- 


-0. 


,0000(0%)=- 


-1 


.5778 


= oE4 


o5= 


0, 


,0469+/- 


-0, 


,047( 


99, 


.87.)+/- 


-0, 


,0000(0%)=- 


-0 


.0133 


= oE5 



From T_Reco_CumObs: : Out_Result_CumObs*** 



cVal Theor MC er_MC(erV/V%) er_st_cor (e/eV%) er_phys (e/eV%) 

A_02/A_P2 1.8000 2.0070 0.072( 3.6%) 0.072( 3.6%) 0.00 ( 0.0%) 

A_P2/A_02 0.5556 0.4983 0.018( 3.6%) 0.018( 3.6%) 0.00 ( 0.0%) 

A_T2/A_02 0.2963 0.3221 0.26 ( 79.5%) 0.012( 3.6%) 0.26 ( 99.9%) 

A02 0.5400 0.5493 0.077( 14.1%) 0.005( 0.9%) 0.077( 99.8%) 

AP2 0.3000 0.2737 0.039( 14.3%) 0.008( 2.9%) 0.039( 97.9%) 

AT2 0.1600 0.1769 0.12 ( 65.5%) 0.006( 3.4%) 0.12 ( 99.9%) 
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cos(d2-dl) -1.0000 -1.0166 0.036( 3.6%) 0.036( 3.6%) 0.00 ( 0.0%) 
cos(dl)*sin(phi) 0.0400 -0.0672 0.14 (201.4%) 0.13 (188.3%) 0.048( 35.5%) 
cos (d2)*sin(phi) -0.0400 0.1350 0.17 (122.7%) 0.13 ( 99.8%) 0.096( 58.2%) 

EUD -2prime- END 

[. . . output for other three observables is placed here] 
SM_DG_TGp_TGO from T_AngMomMethod: : Out .Result _TimeIntObs# BEGIN 



cVal Val 


_TheorC 


Val_MC+/- 


DGL_0 


-0, 


.4566 


-0.4493+/- 


DGL_1 


-0 


.4566 


-0, 


.8170+/- 


DGH_2 


-0, 


.4548 


-0, 


,6876+/- 


DGL_4 


-0, 


.4566 


-0. 


.5358+/- 


DG_0 


-0, 


.4557 


-0, 


.5685+/- 


DG_1 


-0, 


.4557 


-0, 


.7523+/- 


DG_4 


-0 


.4557 


-0, 


.6117+/- 


G_0 


2, 


.2784 


2, 


.2184+/- 


G_l 


2, 


.2784 


2. 


.3103+/- 


G_4 


2, 


.2784 


2. 


.2401+/- 


DG_0/G_0 


-0, 


.2000 


-0. 


.2563+/- 


DG_1/G_1 


-0 


.2000 


-0, 


.3256+/- 


DG_4/G_4 


-0, 


.2000 


-0. 


.2731+/- 



_MC (erV/V%) +/- er_sys (e/eV%) 



0, 


.0388( 


8, 


.6%)+/- 


0, 


.0000 


(0%) 


0, 


. iiisc 


13, 


.6%)+/- 


0, 


.0000 


(0%) 


0, 


.0912( 


13, 


.3%)+/- 


-0, 


.0000 


(0%) 


0, 


.1068( 


19, 


.9%)+/- 


0, 


.0000 


(0%) 


0, 


.0496( 


8, 


.7%)+/- 


0, 


.0000 


(0%) 


0, 


.0720( 


9, 


.6%)+/- 


0, 


.0000 


(0%) 


0, 


.0703( 


11, 


.5%)+/- 


0, 


.0000 


(0%) 


0, 


.0248( 


1, 


.1%)+/- 


0, 


.0000 


(0%) 


0, 


.0360( 


1, 


.6%)+/- 


0, 


.0000 


(0%) 


0, 


.0351( 


1, 


.6%)+/- 


0, 


.0000 


(0%) 


0, 


.0225( 


8, 


.8%)+/- 


-0, 


.0000 


(0%) 


0, 


.0316( 


9, 


.7%)+/- 


-0, 


.0000 


(0%) 


0, 


.0317( 


11, 


.6%)+/- 


-0, 


.0000 


(0%) 



[... output for other two types of DG extraction is placed here] 
END *** T_Accum_Measure: :Average_Measure_OutResult() END *** 
####### T_run_Read: :Loop End T_run_Read: :Loop() . ######## 
***** mainO End ***** 
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Figure captions 

Figure 1. Structure of the package BtoVVana. Inserted squares mean derived classes. 
Diamond inside the class means the data members of the classes shown by arrow. 

Figure 2. Definition of physical angles to describe decays — > J/'0(— > 
l^l^) 0(— K'^K") in the helicity frame [1]. 

Figure 3. Angular distributions obtained by means of class T_run_Read for 
100000 events with S° Jji^ir^ l+l-)(f){-^ K+K') decays ("main" setting of 
physics parameters, see section 5). 
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Figure 2 
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